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Computations are performed to study the boundary layer instability mechanisms 
pertaining to hypersonic flow over blunt capsules. For capsules with ablative heat shields, 
transition may be influenced both by out-gassing associated with surface pyrolysis and the 
resulting modification of surface geometry including the formation of micro-roughness. To 
isolate the effects of out-gassing, this paper examines the stability of canonical boundary 
layer flows over a smooth surface in the presence of gas injection into the boundary layer. 
For a slender cone, the effects of out-gassing on the predominantly second mode instability 
are found to be stabilizing. In contrast, for a blunt capsule flow dominated by first mode 
instability, out-gassing is shown to be destabilizing. Analogous destabilizing effects of out- 
gassing are also noted for both stationary and traveling modes of crossflow instability over a 
blunt sphere-cone configuration at angle of attack. 
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Nomenclature 

distance from cone tip at which mass flux injection is switched on for slender cone case 
boundary-layer similarity length scale (vs/U e ) m 
mass flow rate 
N-factor 

surface distance from the tip or nose of a body along an axisymmetric body generator 

heat flux 

Reynolds number 

streamwise mean flow velocity 

streamwise perturbation velocity 

cross-stream or wall-normal velocity 

streamwise coordinate 

wall normal distance 

wall normal distance normalized with boundary length similarity length scale, / 
density 

angular coordinate 
kinematic viscosity 


Subscripts 

cf - crossflow 

e = boundary layer edge 

n - wall-normal direction 

w = wall 

oo = freestream 
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I. Introduction 

Blunt aeroshells with ablative heat shields are often used to survive the high energy entries into a 
planetary atmosphere. The Orion or Crew Exploration Vehicle (CEV) designed to support human space 
missions returning from low earth orbit and the moon will have to use an ablative heat shield at lunar 
reentry velocities. Whereas stability characteristics and associated transition mechanisms for smooth 
surface boundary layer flows over slender body vehicles are reasonably well understood at this point, 
those for the blunt capsule vehicles are just beginning to be addressed. 1,2 Given its influence on surface 
heat transfer, skin friction drag, and flow separation characteristics, understanding and prediction of 
boundary layer transition over the vehicle surface constitutes an important aspect of hypersonic vehicle 
design. 

A number of experiments have been performed to study the effect of surface ablation on boundary layer 
transition. Typically, most experiments have focused on measuring the effects of surface out- gas sing 3 " 5 
alone or the role of ablation induced distributed roughness, 6 but not both. Roughness induced transition 
in a high-speed boundary layer is a difficult problem. 7 However, an intense fundamental research effort 
funded by NASA and AFRL over the past few years has provided significant clues into the physical 
mechanisms underlying transition due to discrete surface roughness (i.e. see references cited in Ref. 8). 
Extension of these studies to distributed surface roughness is yet to be undertaken. On the computational 
side, aerothermal analysis for ablative heat shields has traditionally been performed without a direct 
coupling between ablation and the material response. However, considerable effort has been devoted in 
recent years towards the eventual development of a coupled capability by implementing surface mass 
transfer boundary conditions into CFD solvers. 9,10 The effects of surface blowing on boundary layer 
instability have also been studied recently in the context of a flat plate boundary layer. 11,12 

This paper is focused on the transition mechanisms due to the effects of out-gassing in the absence of 
surface roughness. Previous studies related to the effects of blowing on flat plate boundary layers 11,12 are 
extended first to a slender cone and then to a blunt capsule. Results for each of these geometries are 
presented in Sections III, IV and V below, following a description of the analysis. For simplicity, the 
flow conditions are such that no high-enthalpy effects are expected in either of the two flow 
configurations. 


II. Analysis Codes 

The unperturbed boundary layer flow over the cone surface was computed on various grids using a 
second order accurate algorithm as implemented in a finite-volume compressible Navier-Stokes flow 
solver VULCAN. 11 The VULCAN computations utilized the code’s built-in capability to accomplish 
shock adaptations. The out-gassing surface boundary condition implemented in VULCAN was based on 
the methodology described in Ref. 10. This methodology requires the specification of the surface normal 
mass flux, gas composition and static temperature. A quadratic equation is solved to obtain the surface 
density and the surface pressure is finally computed via the ideal gas law. The VULCAN implementation 
also included the flexibility to specify an arbitrary out-gassing distribution (or profile) on the boundary 
surface. In selected cases, the mean flow was computed on multiple grids to enable an assessment of grid 
convergence. In select cases, the mean flow was computed on multiple grids to enable an assessment of 
grid convergence. For the case of a slender cone, the mean flow was also computed using a boundary 
layer solver 13 to ensure that the effects of out-gassing on stability characteristics based on the Navier- 
Stokes mean flow compare well with the results based on mean flow obtained by solving boundary layer 
equations. 


* http://vulcan-cfd.larc.nasa.gov 


2 

American Institute of Aeronautics and Astronautics 



The stability of the computed boundary layer flow was analyzed using the Langley Stability and 
Transition Analysis Code (LASTRAC). 14 Analysis was performed using both quasi-parallel stability 
theory and the parabolized stability equations (PSE). Sutherland’s law is assumed to describe the 
viscosity variation for both the mean flow and the unsteady perturbations associated with boundary layer 
instability waves. Stokes law is assumed for bulk viscosity. 

III. Slender Cone 

First, consider the effects of air injection into a Mach 7.4, air boundary layer over a 5-degree semi-angle 
circular cone with a sharp tip. The injection was turned on at s = L, where s denotes the distance from the 
tip along a cone generator and L = 0.09525 meters (3.75 in.) corresponds to the location where the surface 
mass injection begins in the experiment 3 4 modeled herein. The flow conditions and the spatial distribution 
of injection velocity are modeled after Ref. 4. The flow conditions are shown in Table I and perfect gas 
behavior is assumed. 


Table I. Flow conditions for the slender body. 


Mach 

7.4 

Unit Reynolds Number (1/m) 

9.252x1 0 6 

Freestream Temperature (K) 

69.72 

Surface temperature (K) 

308 

Angle of Attack (deg.) 

0 


The computational grids used for most of the axisymmetric mean flow computations (corresponding to 
varying levels of surface mass flux) have 1153 points along the length of the cone and 513 points in the 
wall normal direction. The sharp tip of the cone is approximated by a spherical tip that is 4.572xl0" 6 
meters in radius and is smoothly joined by a 4 th order polynomial to the straight side of the cone. The 
spherical cone tip is resolved by approximately 40 points up to where the sphere tangentially joins the 
straight side. The normal grid is clustered near the wall with approximately 200 points to resolve the 
boundary layer thickness. The length of the cone throughout all computations is fixed at 0.508 meters. A 
denser grid is also constructed with double the number of grid points in each direction to verify grid 
convergence of the mean surface heat flux and boundary layer profiles. Tangential velocity and 
temperature profiles at selected locations along the cone surface are plotted in Fig. 1, which includes the 
results obtained using both of these grids. Virtually no difference can be discerned between the two sets 
of profiles. Mean flow solutions were also obtained using an alternate boundary condition corresponding 
to boundary layer bleed and the profiles (not shown) matched well with the solutions shown in Fig. 1. 

Computations are performed for five selected cases with increasing surface mass flux but an identical 
mass flux distribution along the cone surface based on the measured data as reported in Ref. 4. Case I 
corresponds to an impermeable cone surface (i.e. zero mass injection), whereas cases II through V 
correspond to increasing surface mass flux. The surface mass flux distribution for case II is shown in 
Fig.2. The “amplitude” of the surface mass flux for cases I (zero mass flow), IV and V was selected to 
provide a close match between computed and measured surface heating distribution (Fig. 3). 

Heat flux distributions for five selected magnitudes of the overall surface mass flux are plotted in Fig. 3, 
including a comparison with heat flux data from the experiment. It is seen that, for cases I, IV, and V, the 
surface heat fluxes match those in the experiment. Because of certain ambiguities with the actual mass 
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flux distribution in the experiment, the “amplitude” of the surface mass flux in the computations was 
selected to provide a close match with the surface heating distribution. 

In high Mach number flows, the dominant instability mechanism is associated with a generalized 
inflection point in the boundary layer profile, which is defined as the location where the product of 
density and the normal gradient of streamwise velocity has a local maximum. Fig. 4 shows wall-normal 
profiles of this product at various streamwise locations for cases I, IV and V. The peaks in these profiles 
are clearly seen. For case I, namely, the zero mass flux case, the boundary layer is approximately self- 
similar and, therefore, all profiles nearly collapse onto one another. With finite surface mass flux for cases 
IV and V, the boundary layer thickens and the generalized inflection point moves away from the wall. 

Amplification characteristics of second mode waves for cases I, IV, and V, respectively, are computed for 
a range of frequencies from 100 to 775 kHz in 25 kHz intervals using the linearized version of 
parabolized stability equations (PSE). The logarithmic amplification ratio (N-factor) for axisymmetric 
second mode disturbances at selected frequencies is plotted in Fig. 5 as a function of the streamwise 
coordinate normalized by L, the distance from the cone tip at which the surface mass flux is switched on. 
Results in Fig. 5(a) correspond to the baseline case of an impermeable cone, whereas those in Fig. 5(b) 
and 5(c) correspond to case IV and case V, respectively, from Fig. 3. It is interesting to observe that the 
peak N-factors in Fig. 5(b) and 5(c) are actually lower than that in Fig. 5(a), indicating that the second 
mode disturbances are stabilized by the presence of mass injection in this case. Furthermore, the peak 
frequency is reduced from 250 kHz in case I to 175 kHz in case V, reflecting the effect of boundary layer 
thickening. The physical mechanism behind the stabilizing effect of out-gassing on the second mode 
disturbances has not been identified using current computations. However, one interesting feature of the 
amplification characteristics at nonzero injection corresponds to the extended region of nearly constant N- 
factor at large x (see, for instance, the N-factor curve for the most amplified disturbance of 250 kHz in 
Fig. 5(b), which exhibits an alternating weak decay and weak growth at x/L > 3.5). 

The above results show a stabilizing effect of out-gassing on second mode amplification. The experiment, 
on the other hand, indicated a forward movement in transition with an increasing magnitude of the 
blowing velocity at the surface. This can only mean that the earlier transition in the experiment could not 
have been caused by the second-mode instability mechanism. The question of what exactly led to the 
earlier transition in the experiment remains to be ascertained, especially since it appears unlikely that the 
injection mechanism also acted as a source of unsteady disturbances in the range of second mode 
frequencies, contributing to a significant increase in receptivity for the nonzero injection cases. One other 
possibility that we are investigating at present corresponds to the effect of the porous surface (which is 
not modeled in our current stability equations) on disturbance growth at the frequencies of interest. 


IV. Blunt Capsule 

Findings for a Mach 7.32 flow over a blunt, hemi-spherical capsule are presented next, where, again, we 
assume the injection species to be air, i.e., the same as the incoming flow. The radius of the body is 
0.0889 meters (3.5 in). The flow conditions are listed in Table II. The Mach number at the boundary 
layer edge is subsonic over a significant portion of the body length and the maximum edge Mach numbers 
have low supersonic values. Mach number contours for the baseline case (zero blowing) and the case of 
maximum injection velocity used this study (4 percent blowing) are shown in Figs. 6(a) and 6(b), 
respectively. Besides its obvious impact in substantial thickening of the boundary layer region near the 
surface, the out-gassing also influences the shock location, which gets pushed successively further from 
the body albeit by a small distance. 
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Table II. Flow conditions for a Mach 7.32 hemi-spherical capsule. 


Mach Number 

7.32 

Unit Reynolds Number (1/m) 

14.6xl0 5 6 

Freestream Temperature (K) 

65 

Surface Temperature (K) 

300 

Angle of Attack (deg.) 

0 


The computational accuracy was assessed by using three grids belonging to the same family, but with 
significantly different sizes. An example of the grid convergence test for an injection mass flux of 4% is 
shown in Fig. 7. The coarsest grid has 129 points along the hemi-sphere surface and 353 points in the wall 
normal direction with at least 100 points in the boundary layer. The two finer grids have, respectively, 
twice and four times the grid points in each direction. Comparison of boundary layer profiles based on 
these grids shows that the coarsest grid has slight inaccuracies in this high injection case, while the two 
finer grids yield grid-converged solutions except in the region very close to the outflow boundary (Fig. 7). 
For stability and transition analysis with region of interest far ahead of the outflow boundary, the mid- 
level grid is, therefore, sufficient to ensure accurate computations of the growth of instability waves. 

For the zero mass injection case, the boundary layer flow over the hemispherical capsule does not support 
any significant amplification of instability modes. The absence of instability in the baseline case is 
attributed to the strong favorable pressure gradient over the body surface coupled with a sufficiently low 
Reynolds number. The boundary layer profiles do not have any generalized inflection point, which, as 
explained in Section III, is defined as the location where the product of density and the normal gradient of 
streamwise velocity has a local maximum. Hence, it cannot support any instability modes associated with 
an inviscid mechanism. 

Introduction of mass transfer at the surface produces a dramatic effect on the flow stability when the 
injection rate parameter, m, becomes sufficiently large. Here, m = p w VJpJJ^, where p denotes the fluid 
density, V represents the velocity normal to the wall, subscripts w and oo denote conditions at the wall and 
in free-stream, respectively. At m - 0.01, for example, the previously non-inflectional boundary layer 
profiles have developed a very pronounced peak in the middle of the boundary layer, which moves farther 
away from the surface (both in absolute sense and relative to the boundary layer thickness) at increasing 
distance from the stagnation point (Fig. 8). This modification of mean flow profiles results in the onset of 
inflectional first mode instability. Representative mode shapes associated with the streamwise velocity 
fluctuation produced by this instability are shown in Fig. 9. Consistent with the outward movement in 
inflection point location at increasing distance from the stagnation point, the peak of the (normalized) u 
mode shape moves farther away from the surface as the distance parameter X is increased. Here, u 
denotes the streamwise (i.e., tangential to the surface) velocity fluctuation associated with the instability 
wave. 

The most amplified first mode instability is, in general, an oblique wave. This means that a very large 
parameter range in the frequency and azimuthal wave number space need to be searched to find the 
maximum amplification. However, in this particular case, the boundary-layer edge flow over most of the 
region of interest is subsonic and test runs at a few representative locations on the surface confirm that the 
most amplified first mode waves for all relevant frequencies are very nearly (if not exactly) axisymmetric 
( i.e. having zero azimuthal wavenumber). Therefore, in the first mode instability analysis below, only 
axisymmetric first mode waves are sought, which reduces the parameter range of search by an order of 
magnitude. 
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The effect of out-gassing on the amplification of these inflectional instabilities is shown in Fig. 10, where 
cumulative amplification ratios (N-factors) for selected fixed-frequency disturbances are plotted for two 
values of the out-gassing parameters, m = 0.01 and m = 0.013. All N-factors are computed with the quasi- 
parallel assumption with streamwise surface curvature effects accounted for. Despite the fact that the 
injection velocity increases by only 30 percent from Fig. 10(a) to Fig. 10(b), the associated value of peak 
N-factor is nearly doubled. This suggests that the transition onset location will rapidly shift inward once 
the out-gassing parameter has crossed some threshold value. These observations are qualitatively 
consistent with the measurements by Kaattari. 3 A quantitative comparison with the measured transition 
locations will require using the non-uniform mass flux distribution inferred from the measurements, along 
with an analysis of the impact of the associated uncertainties on the instability wave amplification. These 
computations are currently under way and will be reported in a future paper. 

V. Mars Science Lab (MSL) Capsule 

To examine the effect of surface injection on crossflow instabilities on a blunt configuration, a sphere- 
cone configuration that was tested in the wind tunnel for the Mars Science Laboratory (MSL) is 
considered. 15 The model has a spherical nose of radius 38.1 mm, which smoothly joins a straight 
side of 50 degree semi-cone angle. The base radius of the model is 80 mm. The flow conditions 
are listed in Table III. Again, the perfect gas assumption is used and the mass-flux distribution 
along the front surface is assumed to be uniform except in the immediate vicinity of the shoulder 
region adjacent to the base of the model, where it is tapered to zero. 


Table III. Flow conditions for a Mach 6 MSL capsule model. 


Mach Number 

6.0 

Unit Reynolds Number (1/m) 

15xl0 b 

Freestream Temperature (K) 

60 

Surface Temperature (K) 

300 

Angle of Attack (deg.) 

16 


Two grids belonging to the same family, but of different sizes, are used to verify the grid convergence of 
the mean boundary layer flow over the 50-deg sphere cone configuration. The coarse grid has 97x65x353 
points in the streamwise, azimuthal and wall-normal direction, respectively, whereas the denser grid has 
1.5 times the number of grid points in the streamwise and azimuthal directions, respectively 
(145x97x353). Fig. 11(a) and 11(b) show the boundary layer velocity profiles along the windward and 
leeward symmetry lines. Slight differences are observable at large distance away from the nose of the 
model, but they are not expected to affect the stability analysis within the main region of interest. 

Figs. 12(a) and 12(b) show, respectively, the mean flow density contours on the model wall and the Mach 
contours in the plane of symmetry that contains the windward and leeward lines. The Mach number 
contours indicate the shift in the stagnation point to the lower half of the model. The flow within the 
shock layer region on the lower, windward side corresponds to lower Mach numbers but higher 
temperatures in comparison with the flow on the leeward side. 

One of the important gauges of crossflow instability is the crossflow Reynolds number. The effect of out- 
gassing on crossflow Reynolds number is shown in Fig. 13, demonstrating a significant increase in the 
maximum crossflow Reynolds number as a result of out-gassing, specifically from about 640 to 
approximately 750. 

To further assess the grid convergence for nonzero out-gassing, stability computations were performed 
with the above two grids for the stationary crossflow instability along a streamline that cuts through the 
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region of peak crossflow Reynolds numbers (viz., the 4 th streamline from the leeward symmetry plane in 
Fig. 14 below). The maximum stationary crossflow N factor decreases only slightly from 9.07 for the 
fine grid to 9.01 for the coarse grid, indicating a satisfactory grid convergence for the metric of interest 
herein. It may be noted in passing that a higher grid sensitivity was noted for the amplification 
characteristics of first mode instabilities along the leeward line of symmetry, which is consistent with the 
observations in Ref. 17 for the zero out-gassing case. 

As described in a previous paper, 16 two instability mechanisms dominate in the boundary layer flow over 
this sphere-cone model. First-mode instability is evident along the leeward symmetry plane while 
crossflow instability prevails elsewhere on the leeward side. Fig. 14 shows the effect of out-gassing on 
the evolution of crossflow instability modes along selected streamline paths that were used to integrate 
the amplification rates of stationary and crossflow instabilities along the model surface. The starting 
point of each streamline indicates the onset of instability along that trajectory. The streamlines stop at the 
location where a particular N-factor has been reached. These maximum N-factors for stationary and 
travelling crossflow modes are set to be 5 and 7, respectively. Comparison with the experimentally 
obtained phosphor thermography image is also included in the figure, and it may be seen that the 
observed transition fronts correlate well with these preset N-factor values for the baseline (i.e., zero 
injection) case. The results also highlight the strong effect of out-gassing on the crossflow modes of 
instability. .The transition fronts have shifted substantially farther upstream as a result of the wall- 
blowing. While not discussed herein, the effect of out-gassing on the first mode instabilities along the 
leeward line is strongly destabilizing (even more so in comparison with the crossflow modes shown in 
Fig. 14), analogous to the findings in section IV for the hemi-spherical capsule. It appears that wall 
blowing introduces additional inflectional points on the velocity profiles, which have a direct impact on 
the first-mode type instability. On the other hand, the crossflow instability is mainly due to the 
inflectional instability of the crossflow velocity component, which is not influenced as strongly by the 
out-gassing. 


VI. Conclusions 

Stability computations are carried out for three flow configurations with and without out-gassing in order 
to obtain initial insights concerning the effects of ablation on the stability of the boundary layer flow over 
entry/reentry vehicles. For the second mode dominated axisymmetric boundary layer over a hypersonic 
slender cone, out-gassing actually leads to moderate stabilization of the boundary layer, contrary to the 
measured trend during a previous wind tunnel experiment. This indicates that the transition in the 
experiment was caused by instability mechanism other than second mode. On the other hand, the 
boundary layer flow over the blunt, hemispherical body is dominated by first mode instabilities, which 
become increasingly more unstable as the magnitude of out-gassing velocity is increased, consistent with 
an available set of measurements. A similar destabilizing trend is also noted for a blunt sphere-cone 
model at an angle of incidence which involves crossflow instabilities. Further studies are underway to 
broaden the parametric range of these cases to provide a deeper understanding of the physical 
mechanisms underlying these disparate effects on different type of instability modes, including the 
receptivity and nonlinear stages of the transition process. 
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(a) Tangential velocity profiles. 



Temperature (K) 

(b) temperature profiles. 


Figure 1. Results of grid convergence test for the slender body. 



Figure 2. Mass flux distribution over the length 
of the cone for case II. The mass fluxes for cases 
III, IV and V are, respectively, 2.61, 4.16 and 
9.71 times larger. 



Figure 3. Effect of surface injection on heat 
transfer over a circular cone. Case I 
corresponds to impermeable cone (i.e., zero 
injection), whereas cases II through IV 
correspond to increasing surface mass flux. 
Lines (computations), symbols (data 4 ). 
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(a) Case I (zero out-gassing). All 5 profiles 
approximately collapse on to a single curve due 
to the nearly self-similar nature of the boundary 
layer in the absence of out-gassing. 



(b) Case IV. Out-gassing mass-flow rate: 4.16 
times that of Case II (see Fig.2 above). 



(c) Case V. Out-gassing mass-flow rate: 9.71 times that of Case II (see Fig. 2 above). 


Figure 4. pdvrfdy profiles at selected stations along the slender cone. From top to bottom, cases I, IV 
and V. For each profile, pdxxldy start with some finite value at the wall, reaches at peak value in the 
interior of the flow and decays to zero as y becomes large. Here y is the wall normal distance in self- 
similarity unit (see Nomenclature). 
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(a) Case I from Fig. 3. Peak frequency is 250 (b) Case IV from Fig. 3. With blowing, the 250 

kHz. The 775 kHz curve represents the highest kHz N-factor curve peaks earlier with lower In- 
frequency included in computation. factor than in case I. 


o 

o 

CO 



(c) Case V from Fig. 4. Peak frequency lowers due to thickening of boundary layer. 


Figure 5. Effect of out-gassing on second mode amplification in cone boundary layer. N-factor 
evolution for disturbances of selected frequencies. L is the distance from the stagnation point where 
out-gassing is turned on. The high frequency N-factor curves (e.g. 775 kHz) amplify mostly ahead 
of location where blowing starts and, therefore, remain unaffected by out-gassing. 
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(a) mass flux = 0 (b) mass flux = 4% (b) angular coordinate 

Figure 6. Mach contours and angular coordinate schematic, (a) No out-gassing; (b) 4% out- 
gassing, thickening of boundary layer due to out-gassing is obvious in the picture; (c) graphical 
definition of angular coordinate, 0, referred to in Figs. 7, 8 and 9 below. 
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Figure 7. Results of grid convergence test for the hemi-spherical body. Injection mass flux is 4 % . 
The angular coordinate, 0 , is explained schematically in Fig. 6 (c) above. 



Figure 8. Normalized profile across blunt 
capsule boundary layer at selected streamwise 
stations for m = 0.010. The angular coordinate, 
0 , is explained schematically in Fig. 6 (c) above. 
The y-coordinate is non-dimensional as given in 
Nomenclature. 



Figure 9. Mode shapes of streamwise velocity 
perturbation associated with most amplified 
mode at selected streamwise stations. The 
angular coordinate,#, is explained 
schematically in Fig. 6 (c) above. The y- 
coordinate is non-dimensional as given in 
Nomenclature, m = 0.01. 
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S(m) S (m) 

(a) m = 0.0010 (b) m = 0.0013 


Figure 10. Effect of out-gassing parameter m on instability amplification in boundary layer over 
a blunt capsule. N-factor curves as a function of distance along the surface for disturbances at 
selected frequencies. N-factor computations are carried out under quasi-parallel assumption with 
surface curvature effect accounted for. 



Tangential Velocity (m/s) Tangential Velocity (m/s) 


(a) Windward (b) Leeward 

Figurell. Velocity profiles along the windward and leeward symmetry lines of MSL model. 
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(a) Surface density (b) Mach contours in the plane of symmetry. 

Figure 12. Density and Mach contours of MSL model. 



(a) No out-gassing (b) With out-gassing, m- 0.001 

Figure 13. Effect of out-gassing on crossflow Reynolds number ( R cf ) distribution on the blunt 
sphere-cone model (both cases computed with the fine grid tuned to resolve the respective 

boundary layer flows). 
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(a) No out-gassing (computational) 


Y 



(b) No out-gassing (c) With out-gassing (computational) 
(experimental) m=0.001 


Figure 14. Comparison of predicted and experimental transition fronts using N=5 and N =7 for 
stationary and travelling crossflow instability. The asymmetry shown in (b) is most likely due to 

experimental uncertainties. 
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